%% Plot the probability function for the approximation
pm.rho = 0.738;
pm.theta = 1.15;
l = linspace(2,100,50)';
t = 1 - pm.rho./(pm.theta.*(1-pm.theta+pm.rho)).*l.^(1-pm.theta);

plot(l,t,'k','LineWidth',2)
xlabel('Firm size (L)')
ylabel('t(L)')
saveas(gcf,[filepath.output,'/FigureC2.png'])


       